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Abstract 

In an earlier work[ 1 ] we established that causal Israel-Stewart viscous hydrodynamics is only 
accurate in RHIC applications at very low shear viscosities 4m] s /s < 1.5-2. We show here that 
the region of applicability is significantly reduced if bulk viscosity plays a role in the dynamics. 



1. Introduction 

There has been a lot of recent interest in quantifying the effect of viscosity on observables 
in heavy-ion collisions at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven. Calcula- 
tions are most commonly performed using causal dissipative hydrodynamics|0[3ll3>|5t], though 
co variant transport theory can also be utilizedQHB]- 

Hydrodynamics assumes that the system is near local thermal equilibrium. Its region of 
applicability can only be reliably determined with the help of a fully nonequilibrium theory. In 
a recent work[ 1 ] we used covariant transport to establish the region of validity for the causal 
Israel-Stewart (IS) formulation of viscous hydrodynamics, in a longitudinally boost invariant 
setting with a massless e — 3p equation of state (EOS) and only shear viscosity. We found 
that for typical conditions expected in nuclear collisions at RHIC, IS hydrodynamics is a very 
good approximation (more precisely, 10% accurate in computing dissipative effects) when the 
shear viscosity (77,) to entropy density (.?) ratio is not too large, 4m] s / s < 1.5-2. A useful rule of 
thumb we obtained is that, in order to reach such accuracy, dissipative corrections to pressure and 
entropy must not exceed about 20%. This is only a necessary condition but its main advantage is 
that it can be tested directly from the hydrodynamic calculation. 

Here we study the region of validity of IS hydrodynamics for systems with shear and bulk 
viscosity. Quantifying bulk effects in covariant transport near the hydrodynamic limit is unfor- 
tunately unfeasible, at least for a one-component system with 2 — > 2 scattering, because bulk 
viscosity is at least two orders of magnitude smaller than shear viscosity. In order to proceed, we 
assume that the 20% rule of thumb above applies in the more general "shear+bulk" case as well. 

Because for a massless equation of state bulk viscosity identically vanishes, we here use a 
more realistic result from lattice QCD|8[] (we smoothly merge the EOS onto that of a hadron 
gas at low T < 160 MeV). Reliable calculations of viscosity in QCD are unfortunately not 
available for temperatures relevant for RHIC. We therefore focus on the "minimal viscosity"J2] 
paradigm, i.e., set t] s = s/(47r). For the bulk viscosity, we consider a Lorentzian in temperature 
f = £ m s/[l + (T - T C ) 2 /AT 2 ], where the peak height and width are adjustable parameters. 
Matching to the bulk viscosity estimate from Meyer[10] gives our default parameterization £ m 
0.2. AT = 0.03 GeV and T c = 0.192 GeV. Though our approach is similar to an earlier study by 
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Fries et ai d 1 ill , the main difference is that we also ensure thermodynamic consistency because we 
use the complete set of Israel-Stewart equations of motion. We also map out a much wider range 
of initial conditions. 
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Figure 1 : Entropy produced by the time the system cools down to T = 1 80 MeV (r a 8 fm) relative to the initial entropy, 
plotted as a function of fhermalization time and initial shear stress, for rj s /s = l/(4;r) (left) and 0.5/(4^) (right). 



2. Main results 

We solve the complete set of Israel-Stewart equations (cf. [ 1]) for a boost-invariant Bjorken 
scenario with axial and transverse translational symmetry. We only highlight here the equation 
for bulk pressure 



T dYl _ ( n £) nT n l l d j 
dr \ t J 2 \t di 



(1) 



Here Tn is bulk pressure relaxation time. If the last Urn-.- term is ignored, the entropy generation 
rate per unit rapidity d{dS / drj) / dr = tAt[T1 2 /(£T) + 37rl/(4r] s T)] becomes inconsistent with the 
Israel-Stewart expression for entropy dS /drj = TAjiseq - Tl 2 Tn/{2T() - 37T 2 L T 7r /(4Tj] s )] (n^ is 
the shear correction to the longitudinal pressure, T n is the shear stress relaxation time, At is the 
transverse area of the system, and s eq is the entropy density in local equilibrium). Therefore, 
entropy production was overestimated in hlftl . 

Motivated by kinetic theory we set T n = 6rj s /(sT), but for simplicity take rn = T n (in kinetic 
theory Tn ~ 5t„/3 near the massless limit). Our default initial condition for Au+Au at RHIC 
is eo = 15 GeV/fm 3 (To » 0.297 GeV) at a thermalization time tq = 0.6 fm. For tq = 0.3 and 
1 fm, we scale on an isentropic curve tqsq - const. For initial shear stress and bulk pressure we 
map out wide ranges — po < ttq < po, -po/2 < Hq < po/2 where po is the initial pressure. Three 
choices are of special interest: i) local equilibrium (LTE) no = Ho = 0; ii) Navier-Stokes (NS) 
no = -(4/3)j7,(r )/T , n = -f(7o)/T ; and hi) gluon saturation flj (CGC) n ~ -p , U « 0. 

The main quantity we analyze is the entropy produced relative to the initial entropy AS /So. 
We shall impose AS /So < 0.2 as the condition for region of validity. Because we only consider 
longitudinal expansion, we focus on entropy production until the beginning of the hadronic phase 
Tf = 180 MeV (this is almost identical to entropy production until Tf = 8 fm). Figure [Tfleft) 
shows our results for "minimal" shear and vanishing bulk viscosity. For LTE initial conditions, 
IS hydrodynamics is applicable when to > 0.3 fm. Good accuracy for NS initial conditions 
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requires r > 0.6 fm (n^/po ~ -0.846, -0.585, and -0.488 for t = 0.3, 0.6, and 1 fm), while 
for CGC initial conditions even later tq > 1 fm. 




Figure 2: Relative entropy production until Tr = 180 MeV as a function of initial shear stress and bulk pressure, 
for thermalization times To = 0.3 (top), 0.6 (middle), and 1 fm (bottom). Left column is for f(T) based on Meyer's 
calculation) 1(J], right column is for half that large f(T). Crosses indicate local equilibrium initial conditions, while 
pluses are for Navier-Stokes. In all cases rj s /s = l/(4;r). 

Let us now turn on bulk viscosity. The left column of Figure [2] shows AS /So for our default 
£(T) parameterization. Due to the additional entropy produced, LTE initial conditions now re- 
quire To > 1 fm, while simulations from NS and CGC ones necessitate even later thermalization. 
The situation improves somewhat if bulk viscosity is half as large as our default (( m = 0.1). As 
seen in the right column, LTE initial conditions are then suitable when to > 0.6 fm, while NS 
can be accommodated if to > 1 fm. Results are quite similar if instead the width of the £(T) peak 
is halved, i.e., £ m = 0.2, AT = 0.015 GeV. We checked that shorter relaxation times help only a 

3 



little, as illustrated in Figure|3lleft) for tq = 0.6 fm. On the other hand, the region of applicability 
widens substantially if the shear viscosity is a factor of two smaller (cf. Fig. [fright)). 




Figure 3: Same as Fig.f2]with our default parameters but halved relaxation times (left), or tj s /s = Q.5/(4n) (right). 



3. Conclusions 

We study the region of validity of Israel-Stewart viscous hydrodynamics for conditions ex- 
pected at RHIC, based on the entropy produced during the evolution when both shear and bulk 
viscosity are present. Our results indicate that, unlike it was previously hoped, viscous hy- 
drodynamics does not extend the range of dynamical description to proper times earlier than 
t « 0.6 fm, even for local thermal equilibrium initial conditions - unless the influence of bulk 
viscosity is basically negligible, or the shear viscosity of hot and dense quark-gluon matter is 
significantly below the conjectured "minimal" value of rj s = s/(4n). 
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